Spectral and optical properties of doped graphene with charged impurities in the 

self-consistent Born approximation 
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Spectral and transport properties of doped (or gated) graphene with long range charged impu- 
rities are discussed within the self-consistent Born approximation. It is shown how, for impurity 
concentrations ni mp •> n a finite DOS appears at the Dirac point, the one-particle lifetime no longer 
scales linearly with the Fermi momentum, and the lineshapes in the spectral function become non- 
lorentzian. These behaviors are different from the results calculated within the Born approximation. 
We also calculate the optical conductivity from the Kubo formula by using the self-consistently cal- 
culated spectral function in the presence of charged impurities. 



I. INTRODUCTION 

Graphene, the monolayer allotrope of carbon, has at- 
tracted widespread attention since its isolation^, and 
remains to be the focus of intensive research. Among the 
main reasons for this interest are its remarkable electronic 
properties, described in terms of massless Dirac fermions. 
They make graphene an appealing system to study new 
unconventional physics 3 -, but moreover, easy control of 
the electron density through gating 1 and very high room 
temperature mobilities^ make it also a promising mate- 
rial for applications. Since the first experiments, it was 
acknowledged that understanding the role of disorder was 
essential to describe the electronic properties of graphene 
correctly^—. One of the first issues raised and still un- 
der debate is the origin of the main scattering source 
in graphene transport. The linear-in-density dependence 
of the DC conductivity 2 - was first attributed to charged 
impurities trapped in the substrate^— , which were con- 
sidered the main scattering mechanism. But it was later 
shown that resonant scatterersi^ and ripples^ could also 
account for this linear behavior. Although Coulomb im- 
purities are inevitably present in graphene samples and 
influence the transport properties, they need not always 
be the dominant scattering mechanism limiting the mo- 
bility especially at high densities^. 

The minimal conductivity measured at the Dirac point 
is another experimental observation which can clarify the 
main scattering mechanism in graphene 3 ". While the uni- 
versal value of 4e 2 /7r/i was predicted for uncorrelated 
short range disorder as well as for ripples (described in 
terms of random gauge fields), the experimental value 
was consistently found larger and sample dependent 2 -^. 
These experimental observations can be explained by 
transport theory in the presence of Coulomb disorder and 
electron-hole puddle formatio n 7 i 17 i 18 . DC transport mea- 
surements have therefore provided compelling evidence 
of the relevance of charged impurities in the physics of 
graphene in substrates. 

The comparison of the single particle relaxation time r e 



denning the quantum level broadening and the transport 
scattering time r tr defining the Drude conductivity in 2D 
graphene layers can also be a relevant probe of the na- 
ture of disorder scattering of graphene carriers^. Hong et 
al£Q report, by comparing these two independently mea- 
sured scattering times, that Coulomb impurities play a 
dominant role in graphene samples. But another recent 
experiment of the same type suggests that the main scat- 
tering mechanism in graphene is due to strong (resonant) 
scatterers of a range shorter than the Fermi wavelength 21 . 

The study of the local density of states in samples on a 
substrate 2 — has revealed inhomogeneites which have been 
also attributed to charged impurities 7 although ripples 
may also contribute to this phenomenon^ 3 .. Long range 
charged disorder may also contribute to the broadening 
of the spectral linewidth observed in angle-resolved pho- 
toemission experiments 2 ^ - — . 

The long range charged disorder is also impor- 
tant to understand the recently measured optical 
conductivit y 27 ' 28 . These experiments have revealed de- 
viations from the ideal picture of Dirac fermions which 
have been attributed, at least in part, to the presence of 
Coulomb impurities 2 ^. 

Thus it is expected that Coulomb disorder inevitably 
existing in the environment of graphene will be important 
in many physical properties. In this paper we investigate 
spectral and transport properties of doped (or gated) 
graphene with long range charged impurities within the 
self-consistent Born approximation. 

In general, the Born approximation has been widely 
used to treat the Coulomb disorder—. However, the first 
order Born approximation is restricted by rii mp <C n, 
which is not necessarily satisfied in most experiments. 
More elaborate numerical approache s 7 ' 17 ' 18 ' 31 are not 
hindered by this restriction, but the simplicity of the av- 
eraged theory makes it appealing to extend the known re- 
sults beyond first order. The self-consistent Born approx- 
imation (SCBA) is a non-perturbative approximation 30 , 
and it is the simplest natural way to extend the Born 
approximation. In the case of graphene, this approxi- 
mation has been used mainly for momentum indepen- 
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dent (5-ranged impurities. The more complicated case of 
Coulomb impurities in SCBA has been addressed in the 
evaluation of the DC conductivity onl y 32 ' 33 . In this arti- 
cle, we present the numerical solution of the SCBA equa- 
tions for doped (or gated) graphene, and compute both 
spectral and optical properties by using the calculated 
self-energy and compare our results with previous works 
and with experiments. We observe important changes 
for impurity concentrations rn mp > n, as compared to 
the first order case: a finite DOS at the Dirac point, 
non-linear one-particle lifetimes, non-lorentzian spectral 
functions, and a modified optical conductivity. 

The organization of the paper is as follows: we start 
by reviewing the single Coulomb impurity problem in 
graphene in section [TTJ In section IIIII we discuss the 
treatment of random Coulomb impurities in an aver- 
aged theory, and formulate the SCBA equations. Section 
IIVI describes the results of our work: spectral properties 
and the optical conductivity are discussed. Finally, we 
present our conclusions in section IVl 



II. THE COULOMB IMPURITY PROBLEM IN 
GRAPHENE 

The problem with a single Coulomb impurity has been 
widely studied in both low density limit*— and high 
density limi t 37 ' 41 . This problem is the starting point 
in the description of Coulomb disordered graphene, so 
we highlight its most relevant aspects. The importance 
of the single impurity problem stems from the fact that 
an electron in graphene is affected not only by the bare 
Coulomb potential of the impurity, but also by the other 
electrons that redistribute around it. The effective im- 
purity potential is thus screened by the carriers, and the 
inclusion of the screening is very important to account 
correctly for the effects of a random collection of im- 
purities. Since the exact consideration of the screened 
Coulomb potential is a difficult problem, the screening 
is simplified by assuming a linear response approach. In 
this linear screening limit, the screened Coulomb poten- 
tial is given by 



by 



V(q) 



V (q) 



e(q,U^0) 



(1) 



where e(q, w) is the dynamical dielectric function and can 
be obtained from the density-density correlation function 

as 



l + V (q)U(q,oj) 



(2) 



where Tl{q, uj) is the polarizability of the system. In the 
weak interaction limit the II(q, uj) is calculated within 
random phase approximation (RPA) by summing all bare 
bubbles. 

In ordinary 2DEG the weak interaction limit is defined 
by the interaction parameter r s « 1, where r s is given 



(7m) 1 / 2 ' 



(3) 



where n is the 2D electron density and m is the electron 
effective mass^ 2 - . If the density is lowered, non-linear 
response comes into play™, and the parameter r s deter- 
mines the range of densities where the linear RPA model 
is applicable. 

In the case of graphene, one could expect that the RPA 
is a good approximation at high densities but its valid- 
ity near Dirac point is questionable. In contrast to the 
regular 2DEG, the parameter r s of graphene is density 
independent, and it is known as the coupling constant 
a — e 2 /(KVp). This generates a difficulty in evaluating 
the range of validity of the RPA for graphene. It is un- 
derstood that in the zero doping case it is definitely not 
valid 44 ', while it has been stated that in the doped case 
it is applicable when r s <C 1—, the same criterion for the 
use of normal perturbation theory. Since r s is density 
independent, it is not clear how the two regimes inter- 
polate, and to the best of our knowledge no quantitative 
criterion on n = k F /ir has been established to separate 
them. 

In spite of that, due to its simplicity and its qualita- 
tive success in comparison with experiments, the RPA for 
graphene has been thoroughly studied in the literature 3 . 
In static limit the polarizability is given by— i 



q k F 



(4) 



where qxF = 4a v kp and 
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For k < 2kp, this is just the Thomas-Fermi result 37 , 
and identical to the regular 2DEG 42 . Thus the difference 
in screening between 2DEG and graphene arises at high 
wave vectors, k > 2kp. Even though the high density 
screening is similar to the 2DEG, it is worth noting that 
as n — > RPA screening in graphene only changes the 
constant dielectric constant, but not the wave vector de- 
pendence. This is due to the vanishing density of states 
at the Fermi level as n — > 0. The RPA result is question- 
able at n — 0, where the problem becomes strongly non- 
linear— but its results compare favorably with tight 
binding exact results 4 ^. In this work, both the linear re- 
sponse RPA polarizability and the Thomas-Fermi result 
will be used as a model of screening, bearing in mind the 
previous discussion regarding their limits of applicability. 

Finally, it is also worth noting that when we consider 
many impurities, we face a more complex problem in 
terms of screening (even before disorder averaging), be- 
cause the screening of one impurity may depend on the 
charge accumulated on the rest of them. It is reason- 
able to assume independent screening of impurities when 
the impurity density is small, but this picture may fail 
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when the screening length becomes bigger than the av- where s 



— 1/2 

erage distance between impurities n imp . The many im- 
purity problem introduces another type of non-linearity 
which may become relevant at low densities, where the 
screening is weak 49 . Since the Thomas- Fermi screening 
length is given by 

/TF = ^ = 4a(™)V2' ® 
the condition for independent screening can be written 



as 



i P ^ (4a) 2 7rn. 



(7) 



In the case of graphene on SiC>2, we can take a ~ 0.75 
and this means rii mp /n < 27 which is well satisfied for 
the ranges of n and will consider. 



III. AVERAGED THEORY FOR RANDOM 
CHARGED IMPURITIES 

In the standard theory of disorder— is is assumed that 
the properties of the system with a particular realization 
of the disorder landscape are the same as those averaged 
over all impurity positions. In the case of graphene, the 
low energy Hamiltonian for a particular distribution of 
Ni m p impurities is given by 

H = J d 2 r [vp^adtP - n^iji + ^V(r)ip] , (8) 

where vp = 10 6 m/s is the Fermi velocity, and \x is the 
chemical potential (note that energies are measured with 
respect to the Dirac point). This Hamiltonian is a good 
approximation to the band structure for energies up to 
F> ^5 1 eV—, so we will consider the properties of this 
model in this range of energy. The largest electron den- 
sity we will consider will be n = 6 • 10 12 cm -2 , which 
corresponds to a chemical potential of fi = 0.3 eV, well 
below the limit of applicability. The disorder potential 
V(r) is given by 



N imp 

V(r) =Y1 I d 1 e 



i(r — R,i)q 



V(q), 



(9) 



where Ri are the positions of the Ni mp impurities, and 
V(q) is the Fourier transform of the potential. The disor- 
dered system is described by its Green functions averaged 
over impurity positions, which we assume uncorrelated 
for simplicity. In the case of graphene, this approach has 
been widely used to model disorde r 33 ! 50 " — . When the 
potential is weak, multiple scattering events can be ne- 
glected, and the resulting perturbative series is known as 
the Born approximation, of which only the first terms are 
necessary. The first order term in this series for a general 
potential reads 



E s (fc, w) = 71 



imp 



dk 2 



{-,— for the upper and lower bands, rij mp is 
the density of impurities 

Himp — Ni m pl L , (H) 

with L the system size. Go jS ' are the bare Green functions 
for each band 



G ,±(k,uj) 



u> ± vf \k\ + ie 
and F ss i are the overlap factors 

F SS ' = ^(1 + ss' cos k ,k')- 



(12) 



(13) 



In the self-consistent Born approximation the imagi- 
nary part of the self-energy in the bare Green function 
is replaced by the self energy of the full Green function, 
i.e., ie — > S, and we have the self-consistent equation as 



(14) 



where the self energy is given by 



S s (fc, w) = n 



imp 



dk 2 

W) 



-V 2 {k-k>)J2Gs'(k',u>)F ss ,(k,k'). 

(15) 



The SCBA for the short range case gives well known re- 
sults. Note that, since the potential is independent of k, 
the self-energy in this approximation is also independent 
of k. At u> = it can be solved analytically, giving the 
well known finite purely imaginary self-energy 



S(cj = 0) = 



iA 



exp(: 



and the DOS 



Po 



7m 
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imp v o 



exp(- 



0-i 



(16) 



(17) 



-V 2 {k-k')Y^Go,Ak\uj)F ss ,{k, 



k% 



where A is a high energy cutoff. For undoped graphene, 
however, it has been long known that the SCBA is in 
principle not a good approximation due to the absence 
of a EpT type of parameter that would allow to neglect 
the crossing diagram s 53 ! 54 (but see ref. [55| for a special 
case where the non-crossing approximation is controlled 
by 1 /N, N being the number of valleys) . Renormalization 
group calculations as well as exact results^ indeed differ 
qualitatively from the SCBA. 

Nevertheless, it still represents a good approximation 
in the doped case. While the short range case is ana- 
lytically tractable, the SCBA has not been applied to 
Coulomb impurities because its implementation is not as 
simple 5 ^. The potential due to screened Coulomb impu- 
rities is given by 



(10) 



V C (q) 



2ne 2 1 
Ke(q) q 1 



(18) 
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as stated in the previous section, and the dielectric func- 
tion is given by Eqs. Q and (J5J) (for comparison, we will 
also use the Thomas-Fermi dielectric function, the low q 
limit of Eq. ©). 

In the first order Born approximation, the imaginary 
part of the self-energy has been obtained analytically. 
The computation of Eq. (fit)]) with the Coulomb potential 
Eq. (HHJ) gives, at k=<M: 



which gives 

im[E (1 )(0,cj)] = n lmp ir 2 a 2 



(\k'\+q TF ) 2 2 

(19) 



(u)/v F + q T F) 2 ' 



(20) 



Since this was computed in the Thomas-Fermi approxi- 
mation, it is only valid for uj < 2uip. At k — kp,ui = ujp 
we have 1 ^ 



E = 



a 2 n imp vpTTl(2a) 
kp 



(21) 



where /(2a) is a dimensionless function which can also be 
found in ref. M (for our purposes we will take a = 0.75 
and we have 1(2 * 0.75) = 0.224). This self-energy al- 
lowed to compute the density of states in this approxi- 
mation, which was shown to vanish as |.E|ln|iS| at the 
Dirac points. The SCBA equations for Dirac fermions 
in the presence of Coulomb impurities are obtained by 
substitution of the Coulomb potential Eq. ([T8"l) in Eq. 
(ITS1) . obtaining explicitly 



S+(fc, u) = n imp J [Vc{k - k')f i 
-^-^ [Vc(k-k')} 2 - 



1 + cos 9 1 — cos 9 

uj — vpk' — E+(fc' ,uj) ui + vpk' — £_ (fc', uj) 

1 — cos 9 1 + cos 9 

uj — vpk' — E + (fc', uj) uj + vpk' — £_ (k' , uj) 



Z+(x,y) 



ir n 



x'dx'd9 



\x-x'\ + ^ 



1 + cos 9 



+ 



1 — cos ( 



v F x' -Y, + (x',y) y + v F x' -T,-(x',y) 



(22) 
(23) 



If all quantities with dimensions are scaled with the Fermi momentum or energy, x — k/kp, y = u> /ujf, and S = T,/ujf 
this formula can be rewritten as 

2 



(24) 



with 



kr 



(hit) 



1/2 



(25) 



The iteration of formula Eq. (|24|) generates the series of 
diagrams for the SCBA, and therefore the parameter 
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a 2 rii 



(26) 



controls the relevance of higher order terms. Self consis- 
tent effects become important when 7 ~ 1. Note that 
the independent screening condition Eq. ([7]) in terms 
of this parameter reads 7 < (2a) 4 ~ 5. In real experi- 
ments, doping through a gate allows to reach values of 
n as high as 10 13 cm -2 (see ref. [l|), and the concentra- 
tion of impurities in samples is estimated to reach up to 
nimp = 5 • 10 12 cm -2 in the most disordered one o 16 ' 58 . 

The solution of Eq. (|24p can be obtained numerically 
by discretizing k' and iterating it until an error bound 
is reached. Since the momentum lattice has to be kept 
fixed so the self-energy from one iteration can be fed to 
the next one, a simple Simpson rule for the integration 
proved to be the most efficient way (interpolating in k- 
space at each step allowed for more efficient integration 
algorithms, but the overall performance of this strategy 
turned out to be worse than the simple Simpson rule). 



For the integration in k-space to be reliable, the step in 
the discretization has to be much smaller than the imag- 
inary part of the self-energy, because this determines the 
size over which the function to integrate is significantly 
different from zero. This sets a practical limitation to the 
lowest value of ni mp /n to which we have access, which is 
of the order of 0.1. 

Once the self-energy has been obtained, the spectral 
properties of the system are easily computable. The den- 
sity of states is obtained from 

p{uj) = - lim Im[G(x,x',uj)}. (27) 

The inverse quantum lifetime is defined as T~ 1 (k,ui) — 
/m[E(fc, uj)], and it is in general a function of both w and 
k. Its most interesting value is the on-shell lifetime at the 
Fermi energy T~ 1 (kp 1 Ep). Finally, the spectral function 
is computed as the imaginary part of the Green function 



with 



A+(k,u) 



A(k,uj) = - \A + (k,uj) + A-(k,uj)] , 

7T 



Jm[S + (fc,w)] 



(28) 



(uj - uj+(k) - i?e[E+(fc,o;)]) 2 + 7m[S+(fc, w)] 2 ' 
and similarly for A^(k,ui), where uj±(k) — ±vp\k\. 
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IV. RESULTS 
A. Spectral properties 

We now discuss the self-energies and related spectral 
properties obtained from Eq. (f212j). We start from the 
simple case of the zero momentum self-energy, E(/c = 
0, uj), often taken as an approximation for the self-energy 
around the Dirac point. Fig. [T] shows the imaginary part 
of the zero-momentum self-energy for ni mp — 6 • 10 12 
cm" 2 and several densities. For comparison, the first 



order Born approximation results (i.e., Eq. ([20)) ) are also 
plotted for the same values of the parameters. 
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FIG. 1. (Color online) Imaginary part of the zero-momentum 
self-energy, E(0, w), for n; mp = 6 ■ 10 12 cm -2 and different 
densities, n = 6 • 10 11 , 1.7 ■ 10 12 , 6.7 ■ 10 12 cm" 2 , from top to 
bottom (solid lines). For comparison, the first order analytical 
result given by Eq. (|20[1 is also plotted for the same values of 
the parameters (dashed lines). 
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FIG. 2. (Color online) Different approximations for the self- 
energy, at fixed impurity concentration of 6 • 10 12 cm" 2 , 
and for different electron densities, 6.7 • 10 12 cm" 2 (upper 
panel) and 6 • 10 11 cm" 2 (lower panel). The self-energies are: 
£(fc,e(fc)) (thick line), E(k, k) (dashed line), £(0,u) (dotted 
line). T, TF (k, e(k)), as computed with the Thomas-Fermi ap- 
proximation, is shown for comparison in blue (upper curve). 



Several features should be noted in Fig. [TJ When 
ni m p "C n, the SCBA result agree with the first or- 
der Born approximation result as expected. But when 
the impurity density is comparable with electron den- 
sity, n ~ n,i m p, the disagreement between BA and SCBA 
is measurable, especially at high energies. When ni mp 
is further increased (ni mp > n), the deviation becomes 
significant, and more interestingly a finite value of self 
energy is developed at zero energy, which gives rise to 
a finite density of states at the Dirac point (see below). 
The finite density of states at zero energy is one of the 
main characteristics of the SCBA result as compared to 
the first order Born approximation, which predicts a van- 
ishing DOS at zero energy. 

As discussed in the introduction, the self-energy for 
the short range potentials does not depend on the wave 
vector. Therefore, the simple calculation at k = can 
be used for all wave vectors. However, when the self 
energy is a function of the wave vector and energy as for 
Coulomb disorder, the k = self-energy may be useful to 
analyze the physics near the Dirac point at high densities 
but, in general, it is not a good approximation to analyze 
it at finite wave vectors. Instead of zero wave vector self- 
energy, the most relevant value of the self-energy is at the 



on-shell self energy, E(fc, e(k)), where the Green functions 
are peaked. A first approximation to the on-shell self- 
energy is obtained simply by taking Y,(k,vpk), but the 
true on-shell self-energy should be computed with the 
renormalized dispersion relation, H(k,e(k)) (this is only 
relevant when the dispersion relation is greatly changed 
by disorder). 

The imaginary part of the three different self-energies 
£(fc, e(fc)), E(fc, v F k), E(0, lo) is displayed in Fig.|2]for two 
different impurity concentrations. We find that within 
SCBA the calculated self-energy for all different cases 
shows a finite value at the Dirac point for ni mp > n. 
Aside from this, we observe that the k — self-energy is 
indeed very different from the on-shell self energy. One 
should therefore be cautious when applying this approx- 
imation for general computations. We also observe that 
the true on-shell self-energy E(fc, e(k)) presents almost no 
difference from the approximation E(fc, upfc), which can 
thus be employed safely. Another approximation com- 
monly used in the literature is to employ Thomas-Fermi 
screening for the Coulomb impurities instead of the full 
RPA dielectric function. We also show in Fig. [5]the imag- 
inary part of self-energy E(fc, vpk) with the TF screening 
function. A notable difference is observed compared with 
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the RPA calculation. In particular, note that these two 
curves differ even for uj < 2vpkp. 
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FIG. 3. (Color online) Densities of states as a function of 
energy, for a fixed doping n = 1.6 • 10 12 cm -2 , and varying 
impurity concentration, n; mp = 3-10 12 , 6-10 12 , 10-10 12 cm -2 . 
The linear behavior for the clean case is shown for comparison 
(thin line) Left inset: Variation of the zero energy DOS upon 
simultaneous increase of n and rii mv (see text). Right inset: 
Zero energy DOS as a function of n imp , for fixed doping of 
n = 3 • 10 12 cm -2 . 



The density of states of the system can be calculated 
with the self energy, Eq. (|2"7)) . Fig. [3] shows the densi- 
ties of states as a function of energy for a fixed electron 
density n and several impurity concentrations. A signif- 
icant finite value is observed at the Dirac point for high 
impurity concentrations, in agreement with the previous 
discussion on the self-energy To quantify better where 
the onset of this finite value occurs, in the inset of Fig. [3] 
we show the DOS at the Dirac point as a function of ni mp , 
for n = 3 • 10 12 cm" 2 . A clear threshold is observed for a 
value of 7 ~ 0.4. This finite value of the density of states 
is also similar to the one obtained in the short range 
case. However, there is a significant difference. In the 
short range case, the analytical calculation predicts that 
this finite value is proportional to the cut-off, the only 
scale with dimensions at w = 0. The density of states, 
therefore, remains constant if the rest of the scales in the 
problem are changed simultaneously. To check weather 
this behavior is present in the case of Coulomb impuri- 
ties, we computed the zero energy density of states for 
n-imp — 6n sweeping the value of n. This is shown in 
the left inset of fig. [3l A clear monotonous behavior is 
observed, indicating a more complicated dependence on 
these scales than the constant behavior of the short range 
case. 



Another spectral property of interest is the one par- 
ticle relaxation time r e , given as the inverse imaginary 
part of the on-shell self-energy at the Fermi level. Fig. 
2] shows a plot of the dependence of r e with the Fermi 
level hp and the first order Born approximation results, 
Eq. (|21|) . For each concentration of impurities, the fermi 
momentum \if — a^/ni mp which corresponds to 7 = 1 
and separates the high and low doping region is plotted 
as a horizontal line. It is clearly seen that in the high 
doping region, the SCBA is again indistinguishable from 
the first order result. We also see how in the low dop- 
ing region, deviations from the linear behavior are clearly 
identified. 
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FIG. 4. (Color online) One particle lifetime r e as a function of 
kp, for different impurity concentrations (the physical value 
of h is restored in the units). The full line is the first order 
prediction given by Eq. (|2I[) . The dashed line is a guide to 
the eye. A vertical dotted-dashed line is plotted at the value 
kF = oty/riimp which corresponds to 7 = 1, for each n,i mp . 
Inset: Same plot for larger values of fcp, where the linear 
behavior predicted by Eq. (|21f) is fully appreciated. 

Finally, we discuss the spectral function of the sys- 
tem, given by Eq. (|28|) . The spectral function is a 
physical magnitude that is directly measured in ARPES 
experiment a 24 ' 25 i 59 i 60 . A customary way of representing 
the spectral function is through two sets of plots: its 
constant u> sections as a function of k, known as mo- 
mentum distribution curves (MDCs) and its constant k 
sections as a function of u, known as energy distribution 
curves (EDCs). Here we show both of them in Fig.JSJ for 
n lmp = 3 • 10 12 cm" 2 and n = 2 • 10 11 and 6 • 10 12 cm" 2 . 
We observe the typical broadening of the quasiparticle 
Lorentzian peaks due to disorder, whose width increases 
with increasing impurity density. However, a closer look 
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(a) (b) 

FIG. 5. (Color online) a) EDC at a fixed impurity density of 3 ■ 10 12 cm -2 , for two values of the doping n = 6 • 10 12 , 2 • 10 11 
cm -2 , b) MDCs for the same parameters. Note the change in absolute scale from left to right. 



reveals an unexpected feature. Fig. |6] shows the evolu- 
tion of the MDCs with decreasing doping, and the best 
Lorentzian fit to each curve is also shown. It is appreci- 
ated that as doping is decreased the lineshapes become 
strongly non-Lorentzian. This is due to the momentum 
dependence of the self-energy: the MDCs are obtained 
from Eq. (|29p by fixing the frequency. If the function 
£(£;) changes significantly from its on shell value £(fcrj) 
(with ko defined by ujq = e(ko)) within a scale rj, then re- 
garding the MDCs, the self-energy looks like a constant 
if E(fco,wo) << i], and therefore the MDC looks like a 
Lorentzian. Analyzing our numerical data, it can be seen 
that this condition is fulfilled only for small impurity con- 
centrations, and this is the reason for the non-Lorentzian 
peaks shown in fig. [6] 



A (a.u.) A (a.u.) 




FIG. 6. MDCs at fixed u = 0.8 eV for n imp = 6 ■ 10 12 cm" 2 
and decreasing chemical potential fi = 0.30 eV (top left), 0.20 
eV (top right), 0.09 eV (bottom left), 0.05 eV (bottom right). 
The best Lorentzian fit to the data is shown in red. Note 
how the lineshape becomes manifestly non-Lorentzian as the 
doping is reduced. 



B. The optical conductivity 

In this section we investigate the optical conductiv- 
ity of graphene, which is affected strongly by long-range 
charged impurities. The optical conductivity of the ideal 
(or intrinsic) Dirac fermion model is predicted to be fre- 
quency independent and given by the universal valued 



This frequency independent optical conductivity has 
been measured in both suspended samples™ and samples 
deposited on SiO^. This universal result is observed 
even beyond the energies where the Dirac model is valid, 
due to the smallness of the correction induced by trigo- 
nal warping 6 -^. The absence of sizable electron-electron 
interaction corrections in these experiments has also been 
recently discussed^—. Moreover, this ideal picture may 
be modified by a thermal broadening or a level broaden- 
ing of the single particle states due to disorde r 51 ' 67 . In the 
presence of the broadenings the universal value Eq. (|30[) 
is obtained only in the asymptotically high energy limit, 
where all energy scale of broadenings are negligible. 

Optical conductivity experiments have also been per- 
formed in the gated syste m 27 ' 28 , where the optical con- 
ductivity with different electron densities has been stud- 
ied. Theoretically, the optical conductivity at finite car- 
rier density (or finite chemical potential fi) has been ad- 
dressed by many author o 18 ' 29 ' 56 ' 68 " — . Without any dis- 
order and at T = the optical conductivity at finite 
density arises from two contributions: a Drude peak at 
zero energy from intraband transitions and a constant 
contribution from interband transitions, Eq. (|30l) . start- 
ing at the threshold energy w = 2(1. However, in the 
real experiment 27 , a number of strong deviations from 
these ideal predictions have been observed. plOa One of 
them is a constant background conductivity in the region 
between the Drude peak and the threshold at uj = 2/i, 
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where the interband optical transition is forbidden. An 
anomalously large energy broadening of the threshold 
has also been observed, which cannot be explained in 
terms of thermal broadening. This anomalous optical 
conductivity in the forbidden region has been attributed 
to various physical mechanisms. Short range impurities 
produce a broadening of both the Drude peak and the 
threshold 56 , and more recently Coulomb impurities have 
been considered within Born approximation to produce 
a similar but stronger effect 2 ^. Electron-electron inter- 
actions are also considered to explain the background 
optical conductivit y 71 ' 72 . In addition, the measured op- 
tical conductivity of very low mobility samples in CVD 
graphene 2 ^ shows a reduction of the free carrier Drude 
weight induced by the intraband transition and conse- 
quently a substantial weight increase due to interband 
transition, which can not be explained within Born ap- 
proximation. 

While the effects of short-range scatterers on the op- 
tical conductivity have been computed in the SCBA 56 , 
Coulomb impurities have only been considered within the 
first order Born approximation 2 ^. In this section we cal- 
culate the optical conductivity of graphene within the 
self-consistent Born approximation and provide the ex- 
perimental relevance of our results. We use the Kubo 
formula with the SCBA self-energies calculated in the 
previous section. The optical conductivity can be com- 
puted from the current-current correlation function as 



cr(cj) 



9s9v 
2t:lu 



I dE [ kdk^A s {k 1 E)A s ,(k,E +i 



(31) 

with A s (k,w), the spectral function defined in Eq. (|2U)) . 

Fig. [7] shows the optical conductivity as a function of 
energy for different electron densities with the full self- 
energies obtained from the self-consistent Born approx- 
imation. For comparison, the optical conductivity with 
the first order self-energy given by Eq. (|20[) is also shown. 
As expected, the result approaches to the value <Jq for 
cj » 2fi. Most importantly, the SCBA result gives sig- 
nificantly more background conductivity for the density 
n < riimp. We note that the four curves correspond to 
values of 7 = 0.63, 0.35, 0.22 and 0.16, for given chem- 
ical potentials of fi =0.15, 0.20, 0.25, and 0.30 eV, re- 
spectively. As discussed in the previous section, SCBA 
effects are more significant for the low electron densities. 
The self-energy relevant for this computation is the on- 
shell self-energy, which is clearly different from the k=0 
self-energy. In particular, it gives a non-zero broadening 
at the Dirac point, which is partially responsible for the 
background at u close to fi. 

While these results improve on those computed with 
the first order self-energy for ni mp > n, , they also indi- 
cate that Coulomb impurities are not the only contribu- 
tion in the experiment in^ 7 -. Indeed, any scattering mech- 
anism which is affected by screening will tend to produce 
less scattering as the doping is increased, and what is 
needed to explain the background in the forbbiden region 



is a mechanism which becomes more efficient at high dop- 
ings. Our results are consistent, however, with the latest 
experiment reported imM, which seems to be dominated 
mainly by disorder. 
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FIG. 7. Real part of the optical conductivity for a impurity 



density m- 



6 ■ 10 cm and several chemical potentials 



li = 0.15, 0.20, 0.25, 0.30 eV. The solid (dashed) lines in- 
dicate the optical conductivities calculated with the SCBA 
self-energy (with the first order self-energy at k = 0). The 
corresponding values for 7 for each curve are 0.63, 0.35, 0.22 
and 0.16. 



V. DISCUSSION AND CONCLUSIONS 

We now discuss several aspects of the results shown in 
the previous sections. We have shown that the k = 
self-energy widely employed for the electron self-energy 
is only valid near the Dirac point and very different from 
the on-shell self-energy. Thus it is necessary to consider 
the on-shell self energy to describe the electronic proper- 
ties of disordered graphene. The on-shell self-energy with 
bare single particle energy is not much different from the 
true self-energy and represents a reasonable approxima- 
tion. We also discuss the Thomas-Fermi (TF) screening 
approximation and its relation with the RPA approxima- 
tion. Since the TF screening function is equivalent to the 
RPA result for k < 2kp, it seems that the self-energies 
calculated within both approximations should coincide 
at low energies. However, as shown in Fig. [2] the TF 
self-energy differs from its RPA counterpart at all en- 
ergies. The difference becomes greater at higher values 
of 7 as defined in Eq. (|2"6")) . especially for high impurity 
concentrations. This is because the SCBA is non-local 
in momentum space and the difference between TF and 
RPA at high momenta is enough to alter the low-energy 
region of the self-energy. It is therefore more reliable to 
use the RPA screening function even for uj < 2vpkp. 

We now compare our SCBA results with those ob- 
tained in the same approximation for short range 
disorde r 50 i 51 . We find that even though the densities 
of states for both cases are similar in energy dependence 
at low impurity densities the dependence of both n and 
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n<imp is very different. We also find that the single par- 
ticle relaxation time behaves very differently depending 
on whether 7 > 1 or 7 < 1. The parameter 7 determines 
to what extent the SCBA differs from the first order 
Born approximation. It is a good consistency test that 
for small values of 7, our numerical integration results 
match perfectly with the first order Born approximation 
Eq. (|21l) . This allows us to see clearly the deviation from 
linearity as the disorder density increases. This deviation 
is characteristic of Coulomb impurities which gives rise 
to the ratio r tr /T e to deviate from its constant value at 
low impurity densities. 

The momentum distribution curves (MDC) calculated 
within the self-consistence Born approximation have also 
revealed an unexpected feature. They become non- 
Lorentzian with decreasing electron density, which sug- 
gests that the evolution of MDCs with density could 
be used to probe the relevant scattering mechanisms in 
graphene. If the MDC become non-Lorentzian by low- 
ering the electron density, this is a signature of a mo- 
mentum dependent self-energy, which could be produced 
by Coulomb impurities (other momentum-dependent po- 
tentials such as ripples could perhaps produce a similar 
effect, but the method is still useful to distinguish short 
range scatterers from long range scatterers). 



Finally, the discussed spectral features have a direct in- 
fluence on the optical conductivity. We have shown that 
long range Coulomb impurities beyond the Born approx- 
imation induce a characteristic broadening in the shape 
of the optical conductivity at finite doping compatible 
with recent experiments. 

In summary, this work has addressed the effects of 
long ranged charged impurities on the electronic prop- 
erties of graphene, showing how self-consistent scheme 
modifies them as the impurity density 7ii mp is increased. 
The SCBA in the presence of Coulomb disorder produces 
characteristic features in both the spectral properties and 
the optical conductivity which may be relevant in the ex- 
planation of current experiments. 
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